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1. Introduction 

Consider a rook placed on the lower left corner (0, 0) of a chess board. On how many paths can 
■ the rook reach the upper right corner (n, n) if in a single step it may move an arbitrary number 

of fields upwards or to the right (but not downwards or to the left)? For rectangular chess boards 
of size n x to, the number of paths is given by the coefficient a n m in the rational series expansion 
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The case of square chess boards is consequently the diagonal series of this rational function, which 

O happens to be 
°° 1 1 - x 

XX"*" = 2 + 2V1-10X + 9X2- 

n— 

From here, all sorts of information about the numbers a n ^ n can be easily extracted by means of 
computer algebra, for instance the initial terms 

1, 2, 14, 106, 838, 6802,... (A051708), 

or recurrence equation 

(n + 2)a n+2 ,„+2 - (10n + 14)a n+1 , n+ i + 9n a n ,„ = (n > 0), 

, or the asymptotic formula a n ^ n ~ \jl^3 2rl ~ 1 (n — > oo). Computer algebra can also find the 

algebraic representation of the diagonal series given the bivariate rational series as input, so there 
is altogether no need to do any calculation by hand. 

At a marvelous meeting at Nankai University in August 2010 on the occasion of the second author's 
60th birthday, Frederic Chyzak reported that he and his colleagues had succeeded in doing the 
analogous computation for 3D [3], i.e., they determined the number of paths a rook can take on a 
3D chessboard from (0, 0, 0) to (n, n, n) moving in each step an arbitrary positive integer distance 
into one of the three directions, i.e., moving either by (i, 0, 0) or by (0, i, 0) or by (0, 0, i) for some 
positive integer i. Denoting now the number of this kind of walks by a n , they found the initial 
terms 

1, 6, 222, 9918, 486924, 25267236,... (A144045), 
the recurrence equation 

2(2 + n)(3 + n) 2 (53 + 35n)a„ +3 

- (2 + n)(43362 + 63493n + 30114n 2 + 4655n 3 )a„ +2 
+ (1 + n)(54864 + 100586n + 59889n 2 + 11305n 3 )a„+i 
- 192n 2 (l + n)(88 + 35n)a„ =0 (n> 0), 
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and the asymptotic formula a n ~ 5^64" (n — > oo), all fully rigorous, including so-called cer- 
tificates which allow for an independent formal verification of the obtained results. While it is 
clear in theory that computer algebra is able to obtain this information, it is remarkable that it is 
possible to actually carry out these calculations in practice, because the 3D case requires far more 
computational power than the 2D case. 

If we don't insist on a fully rigorous formal verification, the diagonal recurrence can be obtained 
with much less effort: it suffices to compute some 25 terms of the sequence and use automated 
guessing to find a recurrence which matches them. See [2, 5] for recent developments of this 
technique and the references given there for classical versions. For the present paper we applied 
this technique to empirically find recurrence equations for rook paths in dimensions greater than 
three, and we pose it as a challenge to provide rigorous certificates for them. While at least for 
the very high dimensions this seems totally hopeless for now, we do expect that the coming years 
(or decades?) will see not only faster and bigger computers but also more advanced algorithms 
which can certify our claims within a reasonable amount of computing time. At least we intend to 
encourage progress in this direction. We see no other reason to ask for a certification. The question 
cannot be whether our claims are correct — the empirical evidence is way too strong to leave any 
reasonable doubt about that. Nor can the question be whether there actually exist certificates for 
our claims — it is clear by theory that recurrences of diagonal sequences of multivariate rational 
series can always be certified. Nor can the question be whether a proof may provide some insight or 
understanding — certificates are usually just messy polynomials. The interesting questions instead 
are: how big are the certificates, what is the computational cost for constructing them, and who 
will be the first to get the computation done. 

2. An Alternative Route for Turning our Semi-rigorous proofs to Full-Fledged 

Rigorous Proofs 

We know a priori that there is a recurrence, this follows from general holonomic nonsense. But by 
the work of Moa Apagodu and Zeilberger [1] one can derive a priori upper bounds for the promised 
recurrences. The recurrences for <i-dimensional rook walks turn out (empirically) to have order d 
(for 2 < d < 12). It is very possible that it won't be too hard to prove this sharp upper bound 
in general, or even a weaker — but still realistic — one. This may enable one to give a "soft" proof 
that the empirically "guessed" recurrences are indeed rigorously proved. 

If one would be able to find (realistic!) a priori bounds for the degrees of the coefficients as well, 
then by plain linear algebra the "guessed" recurrences would be rigorously proved. 

3. A short interlude: Why is this problem So important? 

The harsh and/or sceptical critic may say: Who cares? Not that many people (or machines) play 
12-dimcnsional chess, and even the vast majority of the many people that do play traditional 2D, 
8x8 chess, couldn't care less about the number of ways a rook can walk. 

But everyone owes money, and usually to many creditors! The number of ways a rook can walk, 
in the d-dimcnsional cubic lattice, from the origin to [n, . . . , n] is also the number of ways of 
repaying all your creditors if you currently owe n dollars to each of d different creditors, and a 
single payment consists of paying any positive amount of dollars (up to the whole debt) to any 
one of your creditors. Now this is a very practical problem. 

4. Fast Computation of Sufficiently Many Sequence Terms 

As the dimension increases, so does the order of the diagonal recurrence and the degree of the 
polynomials appearing in it. The larger a recurrence is in terms of order and degree, the more 
sequence terms are needed to recover it from sequence data. For instance, in dimension d = 12, 
we needed 1600 diagonal terms in order to find the recurrence. To modern guessing software (we 
used code written by the first author [6]), this is still a moderate problem size. Much harder than 



THE COMPUTATIONAL CHALLENGE OF ENUMERATING HIGH-DIMENSIONAL ROOK WALKS 



3 



guessing the recurrence is the computation of sufficiently many terms on the diagonals, which are 
needed as input for the guesser. The naive way is to start from the rational function 

P(xi,...,x d ) _ 1 - V « r"i T n * 

q[Xi,...,x d ) i 1 _ xi X _ X2 1 _ Xa nu nd=Q 

Its denominator q(x\, . . . ,x d ) gives rise to a multivariate linear recurrence with constant coeffi- 
cients, which can be used to compute the a ni ,..., nd recursively. For example, for d = 2, the rational 
function 

1 = (x-l)(y-l) 

l-I^-I^ l-2x-2y + 3xy 

implies the recurrence 

30n+l,m+l — 2a rim ^i 2a rl _|_i !m + Cl n ,m = 0. 

Together with suitable boundary conditions, this allows the computation of a n , m for arbitrary 
n, m, and hence for a n>n for arbitrary n. 

But this is very costly. In dimension d, in order to compute the nth diagonal term, the recurrence 
forces us to compute all terms a ni ,...,n d with < n\, . . . , n d < n, altogether more than n d terms. If 
n = 1000, a computer won't mind doing this for d = 2, but for d = 3 it is already getting painful, 
and for d > 3 either the memory requirements will exceed the computer's capacity or the runtime 
will exceed the user's patience. Or both. For d > 10 the naive method will not even suffice for 
computing the first n = 10 diagonal terms within a reasonable amount of time. 

Fortunately, there are more efficient recurrence equations. For arbitrary dimension d, we have 

n d a(ni, . . . , n d _ 2 , n d -i,n d ) = {n d -i - l)o(ni, . . . , n d _ 2 , n d _i - 1, n d - 1) 

+ (rid-i + l)o(ni, • • • , n d _ 2 , n d _i + 1, n d - 1) 
+ (2 - n d )a(ni, n d - 2 , n d -i,n d - 2) 

+ (2n d - 2n d _i - 2)a(m, . . . , n d _ 2 , n d -i,n d - 1). 

Note that any application of this recurrence leaves the indices m, . . . ,nd_ 2 fixed, increases n^-i 
and decreases n^. This special form breaks the exponential complexity. It can be shown that 
computing the first n diagonal terms via this recurrence requires only 0(n 2 d 3 ) operations. For 
d < 7, this method was efficient enough to produce enough terms to obtain the recurrence for the 
diagonal. 

For d > 8, an additional improvement was needed. Here instead of directly computing the terms 
on the main diagonal, we first used the previous method for computing the terms of the bivariate 
auxiliary sequence 

up to n, m < 200 or so. Then we used a multivariate guesser to discover some bivariate recurrences 
in n and m for b n _ m and used these to compute the diagonal terms b n _ n = a n .... in . n for n as far as 
needed. 



5. Recurrence Equations 

Most of the recurrences we found are too big to be reproduced here. We make them available 
online at 

http : //www . rise . jku . at/people/mkauers/walks/. 

Here we only give a table with some statistics on their order, the maximal degree of their poly- 
nomial coefficients, and the length of the longest integer appearing in them, measured in decimal 
digits (dd). 
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dim 


ord 


deg 


maxint 


OhlS tag 


comment 


2 


2 


1 


2dd 


A051708 


easy 


3 


3 


4 


6dd 


A144045 


Chyzak et al.'s result 


4 


4 


9 


12 dd 


A181749 




5 


5 


18 


31 dd 


A181750 




6 


6 


31 


51 dd 


A181751 




1-7 

7 


7 


50 


94 dd 


A181752 




8 


8 


75 


149 dd 


A181754 




9 


9 


108 


236 dd 


A181725 




10 


10 


149 


306 dd 


A181726 




11 


11 


200 


462 dd 


A181727 




12 


12 


261 


609 dd 


A181728 





6. Queens 

We have also computed recurrences for the analogous problem of Queen walks, but so far we were 
only able to go up to dimension 5. The relevant output can be found in the above-mentioned 
webpage. 

7. Higher Order Asymptotics 

The leading-term asymptotics for diagonals of rook walks has been derived by Martin Erickson, 
Suren Fernando, and Khang Tran [4] using the powerful analytical method of Robin Pemantle and 
Mark C. Wilson [7]. It turns out to be: 

where ad is given by: 

_ d d + 2 
ad ~ (d + 2) d - 1 (d + l) 2 2 d - 1 ' 
This result matches well with the numbers produced by the recurrences we discovered. For in- 
stance, for d — 12 and n = 500000 we find 

^?;, J 'l' ( , ( ? l)J) '"-'-o°o°o^"--- 

A181728(n) 

And this not all. Thanks to the Maple package AsyRec available from 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/AsyRec 

(see [8]) one can get, very easily, from the recurrences, higher-order asymptotics, using the Birkhoff- 
Trjitznisky method. 

The order-10 asymptotic formulas for the sequences for 2 < d < 9 can be gotten from 
http : //www . math . rutgers . edu/~zeilberg/tokhniot/oRookAsymptotics 

that is based on the input file 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/inRookAsymptotics 

that uses AsyRec and of course, the recurrences obtained by the first-named author's computer. 
A cross check with d = 12 and n = 500000 now yields the very convincing quotient 
^ (n7r) (i- d)/2((rf + 1)d)n(1 + D j + + , . , + n j m) 

Al81728(n) 

= 0.999999999999999999999999999999999999999999999999999999999999999963446 . . . 
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where the □ symbol suppresses some explicit rational numbers which are too lengthy to be repro- 
duced here but which can be also found on the website above. 

And this is still not all. By looking at the output of AsyRec for the sequences for specific d, it 
appears that we have the more refined asymptotic expression for the number of rook-walks from 
[0 d ] to [n d ] for fixed, but arbitrary (symbolic!) d 

^r d {n^- d y 2 {{d + i) d ) n 

(d-l)(d+l)(d 3 + 6d 2 + 18rf+12) 1 



( 



+ 



12d(rf + 2) 3 n 
(d - l)(d + l) 2 (d s + lid 7 + 60d e + 168d 5 - 108d 4 - 564rf 3 - 1632rf 2 - 1584rf - 576) 

288d 3 (d + 2) 6 



We leave the rigorous proof of this as another challenge to the reader. 

8. Fixed n, variable dimension 

Let w n (d) be the number of ways a rook can positively walk from [0 d ] to [n d ]. So far, we fixed d 
and let n vary. But what if we fix n and let d vary? Of course u>o(d) = 1 and wi(d) = d\, Sloane's 
A000142. The sequence W2(d) is of more recent vintage, it is Bob Proctor's sequence A105749. 
But a search on Nov. 19, 2010, did not find W3(d) in Sloane, or elsewhere. 
The Maple package RookWalks available from 

http : / /www . math . rutgers . edu/~zeilberg/tokhniot /RookWalks 
handles these sequences, and the webpage 

http : / / www . math . rutgers . edu/~zeilberg/tokhniot/oRookWalks 

lists the first 150 terms of w n (d) for 1 < n < 4, as well as guessed recurrences and implied 
asymptotics. The asymptotic formulas for the individual n (for 1 < n < 4) lead one to conjecture 
that the leading asymptotics for w n (d) as d — > oo is 



We leave the rigorous proof of this as yet another challenge to the reader. 
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